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First principles density lunctional theory DFT+U calculations and experimental neutron diffrac- 
tion structure analyses were used to determine the low-temperature crystallographic and magnetic 
structure of bixbyite a-Mn2C>3. The energies of various magnetic arrangements, calculated from first 
principles, were fit to a cluster-expansion model using a Bayesian method that overcomes a problem 
of underfitting caused by the limited number of input magnetic configurations. The model was used 
to predict the lowest-energy magnetic states. Experimental determination of magnetic structure 
benefited from optimized sample synthesis, which produced crystallite sizes large enough to yield 
a clear splitting of peaks in the neutron powder diffraction patterns, thereby enabling magnetic- 
structure refinements under the correct orthorhombic symmetry. The refinements employed group 
theory to constrain magnetic models. Computational and experimental analyses independently con- 
verged to similar ground states, with identical antiferromagnetic ordering along a principal magnetic 
axis and secondary ordering along a single orthogonal axis, differing only by a phase factor in the 
modulation patterns. The lowest-energy magnetic states are compromise solutions to frustrated 
antiferromagnetic interactions between certain corner-sharing [MnOe] octahedra. 



Introduction 

Manganese is a multivalent element. Each valence 
state has a characteristic oxide or oxides. Mn ion mag- 
netism and electron correlations of Mn d states make it 
challenging to accurately determine the electronic struc- 
ture of Mn oxides using the density functional theory 
(DFT) approach. Franchini et aL[T] explored various 
manganese oxides using DFT and showed that the elec- 
tronic structure as well as the magnetic ground state de- 
pend sensitively on the details of the DFT calculation, 
such as the exchange-correlation functional used. 
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FIG. 1: Topology of bixbyite phase. Mn in dark blue; O in 
light red; unfilled tetrahedral interstitials indicated in white 



The a-Mii203 phase, with the bixbyite structure, is 
particularly challenging to model because of its compli- 
cated, and not yet completely solved, magnetic structure. 
In the bixbyite structure (Fig. [lj , the Mn 3+ ions are octa- 
hedrally coordinated, while the O ions have 4 Mn neigh- 
bors. The bixbyite structure can be viewed as a close- 
packed lattice of Mn with O ions filling 3/4 of the tetra- 
hedral interstitials in a pattern with IaZ symmetry. Be- 
low about 300 K, a-Mn 2 03 transforms from cubic to an 
orthorhombic structure with Pbca symmetry [2J. Exper- 
imentally, the lattice parameters (a,b,c) saturate at low 
temperatures to approximately a= 9.41 A, b = 9.45 A, 
and c = 9.37 A (Refe.Hgj). Geller[5] rationalized the low- 
temperature orthorhombic distortion of bixbyite as a con- 
sequence of a Jahn-Teller instability of [MnOe] octahedra 
toward elongation along any one of the three Cartesian 
axes. In the cubic bixbyite phase, 24 of 32 Mn atoms ex- 
hibit distorted coordination whereas the remaining 8 Mn 
atoms, which occupy fixed-coordinate high-symmetry po- 
sitions, retain regular coordination environments. The 
orthorhombic phase accommodates Jahn-Teller distor- 
tion of the remaining 8 octahedra (see Fig. J2}-Fig. |3J) . 

The magnetic structure of Mn 2 03 has long been of in- 
terest, but it is not completely solved. Computationally, 
Franchini et al. found a preference for either antiferro- 
magnetic or ferromagnetic ordering in a-Mn203, depend- 
ing on the type of DFT exchange-correlation functional 
used pp. Experimentally, Regulski et al.6\ found evidence 
for various antiferromagnetic ordering transitions within 
the orthorhombic phase, which occur without any appar- 
ent change in symmetry. 

Grant et aL[7j suggested that the magnetic ordering 
of orthorhombic a-Mn 2 03 can be predicted from the 
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FIG. 2: Cubic a-Mn203 viewed along b axis; regular octahe- 
dra shown. 




FIG. 3: Orthorhombic a-M^Og: octahedra shown in Fig. [2] 
have undergone Jahn- Teller distortion. 

cubic Ia3 symmetry because the orthorhombic distor- 
tion is small. They proposed a non-collinear ordering 
model with magnetic moments on the Mn 8(6) sites 
(4(a) and 4(6) in Pbca) aligned with the body diag- 
onals of the pseudo-cubic cell and those on the Mn 
24(d) sites (8(c) in Pbca) directed perpendicular to one 
of the two orthogonal 2i screw axes passing through 
each 24(d) site. However, Regulski et al. demonstrated 
that this model is incompatible with the neutron powder 
diffraction data[5]. They further identified an alterna- 
tive, better-fitting, collinear model which featured anti- 
ferromagnetic ordering on each of the five inequivalent 



Mn sublattices of the orthorhombic structure. However, 
the orthorhombic distortion could not be resolved in the 
diffraction patterns used by Regulski et al and, therefore, 
the atomic positions had to be refined according to the 
high-temperature cubic 7a3 structure. No refinements 
of the nuclear and magnetic structures of the magnetic 
phase under the correct orthorhombic symmetry have 
been reported. 

In this paper, we address deficiencies in both the 
computational and experimental studies of the magnetic 
ground state of a-Mi^O^. Computationally, we use accu- 
rate parameterization of the exchange correlation and on- 
site Hubbard parameters at the density functional theory 
DFT + U level in concurrence with a cluster-expansion 
model to investigate candidate ground state models until 
the correct DFT ground state is established. Experimen- 
tally, we use optimized conditions to synthesize a-Mn203 
with crystallite sizes large enough to yield visible splitting 
of reflection peaks in the neutron diffraction patterns of 
the orthorhombic phase, and thus refine the magnetic or- 
dering within orthorhombic symmetry. Both approaches 
give very similar results for the magnetic ordering, sug- 
gesting that the ground state magnetic structure of a- 
M112O3 is largely solved. 

Computational Methods 

First principles density functional theory (DFT) cal- 
culations, as encoded in the VASP soft ware [9], were 
used to calculate the relaxed configurations investi- 
gated here and their electronic structures. The gener- 
alized gradient approximation (GGA) for the exchange- 
correlation functional was used throughout, within 
the "Perdew-Burke-Ernzcrhof revised for solids" or 
"PBEsol" parameterizationfLO] . 

The PBEsol exchange-correlation functional has been 
found to give excellent results compared with experiment 
for the lattice parameters and bulk moduli of both metals 
and nonmetals[TT]. As is generally true for DFT, how- 
ever, calculated band gaps are too small. This error can 
lead to qualitative errors for narrow bandgap materials 
and for materials with magnetic ions, both which are true 
for a-Mn2 03. We compensated for this error by includ- 
ing onsite Coulomb terms (the "GGA+U" approxima- 
tion). In a previous study of MnC-2 phases [T2], we found 
that the experimental volume and bandgap of /3-Mn02 
could both be reproduced using the rotationally invari- 
ant DFT+U of Liechtenstein et al. [13] with an onsite 
Coulomb potential J = 2.8 eV and onsite exchange po- 
tential U — 1.2 cV for Mn d electrons. Remarkably, these 
same values used for the Mn 4+ ion of /?-Mn02 ion were 
found to be transferable to the Mn 3+ ion of a-Mn2 03, 
giving excellent agreement with experiment, as shown be- 
low. Ex post facto investigations of the effects of varying 
the Mn U and J parameters, or adding U or J parame- 
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ters for oxygen, gave little, if any, improvement. 

Sufficient convergence in total energies and lattice pa- 
rameters was achieved with a plane wave cutoff energy of 
500 eV and a 2x2x2 Monkhorst-Pack grid of k-points. 
A 8x8x8 Monkhorst-Pack grid was used for the density 
of states calculation. The magnetic ordering of each of 
the 32 Mn atoms in the unit cell could be set either "up" 
or "down" as desired. Spin-orbit coupling was neglected. 
Only collinear magnetic structures were computed using 
VASP. 

The aim of the computational work was to find the 
DFT ground state magnetic ordering of a-M^Oa. The 
cluster-expansion concept, as developed for interatomic 
alloys[T21[TS], was used to identify candidate states to ex- 
plore. The formal mathematics of the spin-state problem 
is identical to that of the alloy problem. Each site % is 
given a parameter cr-i, where a\ = 1 for spin up (species A 
in the alloy problem) and ai = — 1 for spin down (species 
B in the alloy problem). 

The total energy U of a configuration {ai} is written 

as 



(i) 



where a are the symmetry distinct geometric clusters, J a 
the effective spin interaction parameter for cluster a, and 
£ a the average spin product over all symmetry equivalent 
occurrences of this cluster [T5]. 

In practice, one calculates individual U for various 
configurations, determines the values of J a according to 
some fitting procedure, and then uses these values to es- 
timate the energies of any configuration, including those 
that are not part of the fit. Errors in the determina- 
tion of J a due to the finite set of included energies lead 
to errors in the predictions. While any method for de- 
termining J a should eventually converge to the same re- 
sult given a sufficiently large number of input configu- 
rations, in this work, we employ the cluster expansion 
method formulated by Cockayne and van de Walle[T7]. 
This method uses DFT results to fit a number of J a pa- 
rameters that is much larger than the number of results. 
The mathematical problem of overfitting the data (that 
is, the non-uniqueness of the solution) is controlled by 
using a Bayesian prior to constrain the magnitudes of 
the J a . This method has additional advantages of ex- 
act fitting of all input energies and self-consistent error 
estimates for all predictions. 

All results were calculated within an identical 80-atom, 
32 Mn cell. The cluster interactions that can be deter- 
mined are limited to those contained within one unit cell. 
The model is valid for predictions of other configurations 
within the same cell, but can not be applied to larger 
supercells; thus if the ground state magnetic state has a 
larger periodicity than the crystallographic one, it will 
be missed. 



There are 2 32 collinear magnetic states for the 32 Mn 
atoms in the a-Mn203 unit cell. Time reversal and crys- 
tallographic symmetry reduce the number of symmetry- 
independent configurations to about 3x 10 8 , which, ac- 
cording to a one-to-one correspondence with the number 
of cluster terms[T3], yield about 3x 10 8 unknown clus- 
ter terms. Solving linear sets of equations with order 
10 8 unknowns is not computationally feasible. To sim- 
plify the problem, we truncated the interaction terms 
at fourth order (time reversal symmetry forbids linear 
or cubic terms in the magnetic interactions). There is 
1 constant term, 73 independent pair cluster term and 
4632 four-body cluster terms in our model. For sets of n 
DFT total energy results, we solved n equations in 4706 
unknowns, using a Bayesian prior to weight the param- 
eters and standard singular value techniques for solving 
underdetermined sets of linear equations. 

The prior that we used for the pair terms was 
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where dij is the distance between the Mn at site i and 
j and do rj 3.3 A= \/2ao/4 is the approximate nearest- 
neighbor Mn-Mn distance and A is an unknown constant. 
The prior for the four-body terms was of similar form, 
with 



. , ,\ . i da ,2 , d Q ,2 , d ,2 , d a ,2 , d ,2 , d a ,2 
W ( hJ ,k,l) = A(-) (-) (-) (— ) (-) (-) 

(4) 

The physical motivation behind the form of this inter- 
action was to weaken cluster terms involving Mn ions 
that are farther apart from each other. Although su- 
perexchange spin interactions are short range, we expect 
that strain coupling effects may mediate longer range in- 
teraction; thus the (dij)~ 2 form for our relative inter- 
action terms in the prior. The value of A was deter- 
mined self-consistently by the leave-one-out cross valida- 
tion method[T7]. This value was scaled such that the root 
mean square error in the predicted energies equalled the 
root mean square of the predicted errors [17 . 

We studied the magnetic states in an iterative man- 
ner. After calculating an initial set of energies versus 
magnetic orderings for a few simple configurations, ad- 
ditional structures were investigated, with, in rotating 
turns, (1) the minimum predicted energy among untested 
structures, (2) the maximum predicted energy among 
untested structures, and (3) the maximum predicted un- 
certainty in energy. The parameters J a were recalcu- 
lated after each step, and then used to predict the ener- 
gies and energy uncertainties for all order 3x 10 sym- 
metry independent configurations. The model was re- 
fined iteratively until there were no more predicted states 
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within two standard deviations of uncertainty of the 
tenth lowest-energy state found, at which point it was 
concluded that the collinear ground state was probably 
found. 76 structures in all were calculated. 

Experimental Methods 

The a-Mn2 03 powder sample was prepared by heat- 
ing MnCC>3 (analytical reagent) at 800 °C in air for 12 
h, which was the highest temperature to yield phase- 
pure a-Mn203 devoid of MJI3O4 traces. The sample was 
characterized using X-ray powder diffraction in an in- 
strument equipped with an incident-beam monochroma- 
tor (Cu K Q i radiation) and a position sensitive detec- 
tor. The heating temperature and time were selected to 
minimize the width of the 222 peak, which remains non- 
split in the orthorhombic phase and, therefore, reflects 
the size of the coherently scattering domains in the sam- 
ple. (Phase-pure a-Mn 2 03 can be obtained by heating 
MnC0 3 at temperatures between 600 °C and 800 °C, 
but lower temperatures produced considerably broader 
peaks). No changes in the peak widths were observed 
after the second heating at 800 °C for 12 h. 

Neutron powder diffraction measurements were per- 
formed using both the time-of-fiight HIPD diffractome- 
ter at the Lujan Center of the Los Alamos National 
Laboratory and the BT-1 constant-wavelength (Cu 311 
monochromator, A = 1.5405 A, 15' collimation) diffrac- 
tometer at the NIST National Center for Neutron Re- 
search. For these measurements, the a-Mn203 powder 
was loaded in vanadium cans. In each experiment the 
data were collected at a series of temperatures (HIPD: 
300 K, 200 K, 150 K, 100 K, 60 K, 40 K, 5 K and BT-1: 
300 K, 100 K, 40 K, 10 K, 2 K). Rietveld refinements 
of the nuclear and magnetic structures were performed 
using GSASHS|. 

The magnetic-structure models were selected ac- 
cording to representational analyses performed by 
SARAhpJ)]; likewise, SARAh was used for symmetry- 
constrained refinements in GSAS. First, the magnetic 
basis-vector coefficients were refined using a Reverse 
Monte Carlo (RMC) algorithm implemented in SARAh 
with the magnitudes of all the magnetic moments con- 
strained to be equal. The best-fit model was further re- 
fined in GSAS (i.e. using least squares minimization) by 
keeping the basis-vector coefficients fixed but allowing 
for distinct ordered magnetic moments on inequivalent 
Mn sites. The HIPD and BT-1 data produced consistent 
structural parameters. 

Computational Results 

The a-Mn2 03 bixbyite structure was first investigated 
with ferromagnetic ordering. Relaxation under cubic 



TABLE I: DFT coordinates for ferromagnetic Q-Mn203 in 
cubic Iai phase. Lattice constant ao = 9.4090 A. 

Species Site x y z 
Mn(l) 8(a) 
Mn(2) 24(6) 1/4 0.2848 
0(1) 24(b) 0.4162 0.1286 0.3555 
Q(2) 24(b) 0.3714 0.1445 0.0838 



TABLE II: Calculated DFT crystal structures for orthorhom- 
bic ferromagnetic and ground state antiferromagnetic struc- 
tures. The similarity of a and b for the orthorhombic ferro- 
magnetic structure is coincidental. 
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Ia3 symmetry, yielded a — 9.409 A, and the struc- 
ture shown in Table [T] DFT phonon results of this cu- 
bic structure show an extremely strong double instabil- 
ity (v = 510 i cm" 1 ), associated with Jahn- Teller dis- 
tortions of the oxygen octahedra centered on the Mn(l) 
sites. Full relaxation of the bixbyite structure perturbed 
by either mode in the instability doublet, or any linear 
combination of the two, leads to an orthorhombic min- 
imum energy state with Pbca symmetry, explaining the 
experimental cubic-orthorhombic transition. In fact, we 
find that all such combinations relax to the same fer- 
romagnetic ground state structure, differing only with 
respect to (1) which of the original cubic axis becomes 
the short axis of the orthorhombic unit cell and (2) pos- 
sible translations of the origin by (1/2,1/2,1/2). Each 
Mn(l) and Mn(2) in the orthorhombic structure has four 
short Mn-0 distances and two long Mn-0 distances. The 
topology of the orthorhombic structure is uniquely de- 
fined by specifying which Mn(l)-0 and Mn(2)-0 dis- 
tances are long. In this work, we arbitrarily choose the 
setting where the Mn(l) at (0,0,0) has its far O neighbors 
at approximately ±(0.139, 0.149, -0.093), and the Mn(2) 
at (1/2,1/2,1/2) has its far O neighbors at approximately 
±(0.656,0.412,0.641). The ground state ferromagnetic 
structure is shown in the left-hand side of Table |ll] 
The orthorhombic crystallographic structure of the fer- 



TABLE III: Ground state DFT magnetic state and lowest en- 
ergy noncollinear magnetic state found in a Heisenberg model 
based on pair interactions determined by fits to DFT results. 
Units are fi b ■ Only the relative directions of the magnetic mo- 
ments are determined; my is the magnetism along the only 
axis for the collinear DFT ground state and principal axis 
of the noncollinear Heisenberg model ground state; m± is 
the magnetism along a second axis of the Heisenberg model 
ground state. Mn positions as in Table |TT| spin moments for 
other Mn are related by applying the factors shown in Ta- 
ble ITVl 
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TABLE IV: Effects of Pbca symmetry generators on the mag- 
netization components of the theoretical a-Mn203 magnetic 
structures listed in Table ITTT1 

Generator mu m± 

(x,l/2-y,l/2+z) -1 -1 

(l/2+x,y,l/2-z) -1 -1 

(l/2-x,l/2+y,z) +1 -1 



romagnetic phase was used as the starting part for re- 
laxation of each collinear spin combination {cr^} studied. 
The state that was ultimately picked as the ground state 
was the 36th studied. In total, 76 states were investigated 
before the termination criterion was reached: no new 
structures with predicted energies within two standard 
deviations of the tenth-lowest-energy structure found. 

The calculated ground state collinear magnetic struc- 
ture is shown in Fig. [4] and listed in Table |III| The spin 
moments are obtained by taking the difference in the 
number of spin up and spin down electrons, integrated 
within spheres of radius 1.24 A centered on each Mn, as 
calculated using VASP. The magnitudes are reasonable 
for high-spin Mn 3+ ions, but are not precisely compara- 
ble with experiment because of the artificial partition of 
the cell into spherical volumes. 

We next investigated possible noncollinear magnetism 
using the parameters found in the fit to the computa- 
tional results. The obvious way to extend the model is 
to use the same parameters for the collinear approach, 
but to vectorize the spins in a Heisenberg model ap- 
proach. That is, instead of using Ising-like interaction 
terms Jija i <Tj 1 one uses Heisenberg-like interaction of 
spin "vectors" according to the principle axes of "up" 
spin: JijUi ■ (jj, with exactly the same set of Jij. For 
simplicity, we only included pair interactions in this ap- 
proach. Simulated annealing was used to find the low- 
est energy state of the Heisenberg model. Random ini- 
tial spin configurations always converged to an equivalent 



FIG. 4: Lowest-energy collinear magnetic structure found 
computationally; also collinear model that gives best fit 
(x 2 = 2.54) to experimental results at 2 K. Only Mn atoms 
shown; dark spheres represent spin "up" and light spheres 
spin "down". 
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noncollinear ground state (right side of Table III I , with a 



FIG. 5: Calculated density of states (DOS) for a-Mn203 in 
ground state collinear magnetic structure. 



dominant spin axis direction my , and all secondary spin 
components along the same orthogonal axis m±. 

The band structure for the lowest-energy collinear 
magnetic state found is shown in Fig. [5] The structure 
has some similarities to that of Franchini et aLpQ, us- 
ing the PBE0 approximation, including a band gap at 
the Fermi level (0.6 eV in our case). Our results differ 
from theirs in that our strict antiferromagnetism leads to 
a band structure without a distinction between major- 
ity and minority spins. Furthermore, we find a remark- 
able signature of the combination of octahedral crystal 
field splitting and effects of Jahn- Teller distortion in the 
nature of the highest occupied levels: they are split off 
from the other occupied 3d states. The highest valence 
band is an isolated band associated with 3d electrons of 
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FIG. 6: A trace of the 622 cubic reflection at three different 
temperatures. This reflection appears as a single peak at 300 
K but exhibits pronounced splitting at sub-ambient temper- 
atures. Similar trends are observed for other reflections. 



the Mn(l) and Mn(2). A second split band at around 
E F - 0.8 eV is associated with Mn(3), Mn(4), and Mn(5) 
3d electrons. 



Experimental Results 

Room-temperature diffraction patterns (Fig.RjJl exhibit 
no clear reflection splitting indicative of an orthorhom- 
bic distortion. However, the peaks are broad and the 
cubic IaZ model fits poorly yielding abnormally large 
atomic displacements parameters Ui SO for the oxygen 
atoms, which suggests that the structure is distorted. 
The data (not shown) can be fitted satisfactorily using 
the orthorhombic Pbca structure reported in the litera- 
ture with sensible Ui SO values. The refined lattice distor- 
tion at 300 K differed considerably between the HIPD 
data (b/c = 1.0007) and BT-1 data (b/c = 1.0036), pre- 
sumably because of the close proximity of the phase tran- 
sition which leads to relatively large changes in the lattice 
distortion even for small temperature differences. The 
distortion increases rapidly on cooling to b/c = 1.0085 
at 100 K and the reflection splitting becomes evident 
(Fig. [6J). For these temperatures, the lattice parame- 
ters refined using the HIPD and BT-1 data were in good 
agreement. Temperature dependencies of the lattice pa- 
rameters, unit cell volume, and b/c data are summarized 
in Fig. [7J Fig. [8] displays the experimental and calcu- 
lated diffraction profiles for T=100 K, while Table [V|sum- 
marizes the results of the nuclear-structure refinements 
at 100 K and 2 K. The MnO*6 octahedra exhibit strong 
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Jahn- Teller distortions already at 300 K (Table VI I 



Below 100 K, a series of strong reflections appears at 
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FIG. 7: Temperature dependence of the (a) orthorhombic lat- 
tice parameters and (b) unit-cell volume and b/c ratio which 
characterizes the magnitude of the orthorhombic distortion. 
The behavior of the orthorhombic distortion changes across 
the magnetic transition at approximately 80 K. The error bars 
are within the size of the symbols. 



lower angles, which signifies magnetic ordering (Fig. |9| . 
The b/c ratio decreases slightly below the magnetic tran- 
sition (Fig. [^b)) to 1.0079 at 2 K. The patterns remain 
qualitatively unchanged from 60 K down to 2 K. The 
magnetic reflections can be accounted for by a propa- 
gation vector k=0. All eight symmetry elements of the 
space group Pbca leave this propagation vector invariant. 
Group-theory analysis yields eight irreducible represen- 
tations (IR) (Ti, i=l, 8), each having an order of one. 
Only four of these IRs (Tj, T 3 , T 5 , and r 7 ), are common 
to all five of the inequivalent Mn sites. 

Regulski et aZ.[6 presented their collinear magnetic- 
ordering model in the form of schematic drawings for each 
Mn sublattice; no symmetry analysis was performed. Ac- 
cording to our representational analysis, the ordering 
types for their sublattices 1 (combined Mn 4(a) and 4(b) 
sites), 2, and 4 belong to the T 3 representation. However, 
the ordering on their sublattice 3 is incompatible with 
any of the IRs; that is, the model, at least as presented, is 
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TABLE V: Parameters of the nuclear structures of paramagnetic (T=100 K) and antiferromagnetic (T=2 K) a-Mn2 03 obtained 
by Rietveld refinements using neutron powder diffraction data (BT-f). In both cases, the space group is Pbca (#61). The 
model assumed isotropic atomic displacement parameters (Ui ao ), which were constrained according to the atom type (i.e. Mn 
or O). The refined values were C/ iso (Mn) = 0.0031(2) A 2 and U iso (0) = 0.0049(1) A 2 at 100 K and U iso {Mn) = 0.0021(3) A 2 
and Ui ao (0) = 0.0040(1) A 2 (O) at 2 K. Numbers in parentheses refer to one standard deviation as calculated in GSAS. 







100 K (x 2 = 1.16) 






2 K ( x 2 = 2.54) 




a 




9.4104(1) 






9.4078(1) 




b 




9.4509(1) 






9.4488(1) 




c 




9.3706(1) 






9.3739(1) 




Atom 


Site x 


y 


z 


X 


y 


z 


Mn(l) 4(a) 

















Mn(2) 4(b) 1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


Mn(3) 


8(c) 0.2598(5) 


0.2854(3) 


-0.0105(5) 


0.268(1) 


0.280(1) 


-0.011(1) 


Mn(4) 


8(c) 0.2867(3) 


0.0017(8) 


0.2436(4) 


0.2880(9) 


0.009(2) 


0.246(1) 


Mn(5) 


8(c) 0.0127(5) 


0.2456(6) 


0.2827(3) 


0.012(1) 


0.237(1) 


0.2833(9) 


0(1) 


8(c) 0.4253(3) 


0.1241(4) 


0.3511(4) 


0.4276(8) 


0.1240(9) 


0.3502(8) 


0(2) 


8(c) 0.1399(3) 


0.3509(4) 


0.4068(2) 


0.1395(8) 


0.3513(9) 


0.4074(7) 


0(3) 


8(c) 0.3579(3) 


0.4175(4) 


0.1218(3) 


0.3586(9) 


0.4169(9) 


0.1209(7) 


0(4) 


8(c) 0.0838(3) 


0.3725(4) 


0.1402(3) 


0.0738(7) 


0.3712(9) 


0.1385(9) 


0(5) 


8(c) 0.3797(3) 


0.1477(4) 


0.0789(3) 


0.3821(8) 


0.1509(9) 


0.0803(7) 


0(6) 


8(c) 0.1554(3) 


0.0882(3) 


0.3609(4) 


0.1563(8) 


0.0882(9) 


0.3586(8) 



TABLE VI: Mn-O distances at 300 K and 100 K (in A). 
Numbers in parentheses refer to one standard deviation as 
calculated in GSAS. 





300 K 




100 K 




Atom 










Mn(l) 


2.03(2) (x2) 




1.949(3) (x2) 






2.02(1) (x2) 




1.944(3) (x2) 






1.97(2) (x2) 




2.129(3) (x2) 




Mn(2) 


2.04(2) (x2) 




1.953(3) (x2) 






1.90(1) (x2) 




2.117(3) (x2) 






2.04(2) (x2) 




1.923(3) (x2) 




Mn(3) 


2.17(2) 


1.88(1) 


2.200(5) 


1.879(5) 




2.03(2) 


2.31(2) 


1.987(5) 


2.327(5) 




1.92(1) 


2.00(2) 


1.915(5) 


1.961(5) 


Mn(4) 


2.00(2) 


2.20(2) 


2.033(6) 


2.181(7) 




1.90(1) 


1.96(2) 


1.934(5) 


1.959(6) 




2.35(2) 


1.88(1) 


2.269(7) 


1.859(5) 


Mn(5) 


1.90(2) 


1.96(2) 


1.889(5) 


1.943(5) 




2.26(2) 


1.91(2) 


2.358(6) 


1.916(5) 




1.98(2) 


2.18(2) 


2.026(5) 


2.134(5) 



not compatible with orthorhombic crystallographic sym- 
metry (possibly there is a drawing error). The spin ar- 
rangement for their sublattice 3 can be made compatible 
with the structural symmetry by flipping the directions 
of two spins. The resulting ordering patterns belong to 
either the T3 or Tq representations, depending on which 
two spins are flipped. We fitted both models (i.e. T 3 
and r3 + Tq) to our data. A collinear pattern with mag- 
netic moments directed along one of the orthorhombic 
axes was assumed. The model with all the moments 
parallel to the c-axis produced a superior fit of quality 
comparable to that reported by Regulski et al. How- 
ever, examination of the misfit between the calculated 
and experimental profiles reveals significant discrepan- 



cies for several magnetic reflections that become split in 
the orthorhombic structure (Fig. [9^a)), which indicates 
that the intensity distribution among the split-peak com- 
ponents is incorrect. Conceivably, these deficiencies were 
obscured by the insufficient resolution in the data used 
by Regulski et al. 

As the literature model failed to describe the data, we 
considered collinear models corresponding to other rep- 
resentations (i.e. r*i, Y§ and TV). The Ti model with 
the magnetic moments directed along the a-axis pro- 
vided a satisfactory fit to the neutron data of quality 
far superior to that obtained for any of the T3 mod- 
els (Fig. [9jb)); the collinear Ti models with magnetic 
moments directed along the b and c-axes yielded consid- 
erably worse agreement factors. Models generated ac- 
cording to the r 5 and representations generated poor 
fits and were discarded. Refinements of the magnetic- 
moment magnitudes (m) independently for each Mn site 
significantly improved the fit. The resulting m-values, 
which at 2 K range from approximately 3 /is to 4 hb, 
are consistent with those expected for Mn 3+ ions. The 
collinear model most consistent with the experimental 
results is shown in Table IVHI It is identical to that de- 
termined computationally (Fig. [4]). 

We explored potential deviations from the collinearity 
using the algorithms implemented in SARAh. A rela- 
tively large number of RMC cycles (3,000 for 2 basis- 
vector mixing coefficients per site and 10,000 for 3 coef- 
ficients per site) were found necessary to locate a model 
that provides an adequate fit (a goodness-of-fit x 2 < 3) to 
the data. Detailed refinements of the magnetic structure 
were performed at 40 K and 2 K. Multiple refinements 
that start from randomly selected values of the mixing 
coefficients were run to verify the robustness of the best- 
fit structural model at 2 K. 
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FIG. 8: Experimental (red/dots) and calculated (blue/line) 
neutron diffraction profiles (BT-1) for a-Mn203 at 100 K. 
The residual is indicated below (green line). The agreement 
factors are \ 2 — 1-16 and R wp = 4.79%. 



The best fit was obtained by restricting the magnetic 
moments to the orthorhombic (010) plane. The fit of the 
several magnetic reflections was visibly improved in the 
non-collinear model (Fig. pic), Fig. 10). The best-fit Fi 
models at both 40 K and 2 K feature similar patterns of 
magnetic ordering ( Table fv|). No additional improvement 
was obtained by varying all three mixing coefficients per 
site after 20,000 RMC cycle; possibly, even this number 
of cycles was insufficient to identify a global minimum. In 
the best-fit non-collinear model, the magnetic moments 
are aligned preferentially with the a-axis but exhibit sig- 
nificant (up to 32°) deviations from this direction; the de- 
viations from collinearity become especially pronounced 
at 2 K. 
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FIG. 9: A low-angle portion of the neutron diffraction pattern 
collected at 2 K showing experimental (red/dots) and calcu- 
lated (blue/line) neutron diffraction profiles. Note that all the 
reflections with 26 < 32° were absent at 100 K (Fig. [If; these 
reflections, all indexable according to the primitive nuclear- 
structure unit cell, originate from magnetic ordering. The cal- 
culated profiles correspond to (top) the F3 model by Regulski 
et al. with magnetic moments aligned with the c-axis, (mid- 
dle) the Ti model with magnetic moments collinear with the 
a-axis, and (bottom) the Fi model with non-collinear mag- 
netic moments residing in the ac plane. The Fi models pro- 
vide a superior fit relative to the T3 model. A non-collinear 
alignment of magnetic moments significantly improves the fit 
for several magnetic reflections (outlines using a dashed line) . 
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TABLE VII: Magnetic-moment (m) components for the Mn 
atoms in the best-fit Ti models at 2 K and 40 K. The atomic 
coordinates at 2 K are given in Table [V] The magnetic- 
ordering model assumed m y =0. The refined magnitudes of 
m (units of fis) are indicated for the inequivalent Mn sites. 
Magnetic moments for the remaining Mn sites are generated 
as described in Table I Villi The refinements at 2 K and 40 
K were performed independently with the random starting 
models for the basis-vector mixing coefficients. The angle be- 
tween the mtot and the a-axis is defined as 4> (°)- Numbers in 
parentheses refer to one standard deviation as calculated in 
GSAS. 



Atom 
Mn(l) 
Mn(2) 
Mn(3) 
Mn(4) 
Mn(5) 



2 K 

m z mtot 



-2.6 
-3.4 
3.2 
3.0 
-3.5 



1.6 3.1(1) 32 

-0.7 3.5(1) 12 

-1.4 3.5(1) 23 

1.3 3.3(1) 24 

-2.3 4.2(1) 34 



-2.3 
-3.0 
3.0 
2.9 
-3.4 



40 K 
in z 
1.4 
-0.8 
-0.3 
0.3 
-1.0 



mtot <p 

2.7(1) 32 

3.1(1) 15 

3.0(1) 5 

2.9(1) 6 

3.5(1) 16 



TABLE VIII: Effects of generators of Pbca symmetry on 
magnetic moments refined experimentally. 



Generator 
(x,l/2-y,l/2+z) 
(l/2+x,y,l/2-z) 
(l/2-x,l/2+y,z) 



-1 
+1 
-1 



Discussion 

The theoretical lattice parameters of orthorhombic a- 
M112O3 change significantly when the magnetic struc- 
ture changes from ferromagnetic to the antiferromagnetic 
lowest-energy collinear state (Table demonstrating 
significant spin-strain coupling. The ground state lat- 
tice parameters are within about 0.1% of experiment, 
phenomenally good agreement that demonstrates the ac- 
curacy of DFT calculations using the PBEsol exchange 
correlational along with on-site U and J parameters to 
treat e?-electron correlations in Mn. 

The lowest DFT collinear state found and the best 
low-temperature experimental fit to a collinear model are 
identical, although they were achieved completely inde- 
pendently, suggesting that the nature of the magnetism 
in a-Mn203 is largely solved. The secondary components 
of the best noncollinear spin arrangements of model and 



experiment (Table III and Table VII I appear different at 



first glance, but in fact, the two are related by the approx- 
imate relationship m± (model) cos(27rx)TO_L(expt.), 
with x the crystallographic positional coordinate. The 
source of the discrepancy is not clear, but the agreement 
is noteworthy given the approximations involved in the 
computational approach. Note that the experimental re- 
sults show clear preferences for magnetic moments along 
particular directions; the DFT approach neglected spin- 
orbit coupling, and thus the effect of magnetic moment 
direction could not be studied. 
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FIG. 10: Experimental (red/dots) and calculated (blue/line) 
neutron diffraction profiles (BT-1) for Q-Mn203 at 40 K (top) 
and 2 K (bottom). The calculated profiles correspond to the 
Fi model with the non-collinear array of magnetic moments 
in the ac plane. The agreement factors are \ 2 = 1-22 and 
Rujp = 5.21% (40 K) and x 2 = 2.53 and R wp = 6.89% (2 K). 
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FIG. 11: Calculated magnetic interaction parameters between 
Mn ions in a-Mn2 03. Positive values for J favor antifer- 
romagnetic alignment. Mn subscripts refer to the different 
Mn sites in the orthorhombic phase. Error bars indicate plus 
and minus one standard deviation of the parameter, based on 
cross-validation calculations. 



The pair magnetic parameters, as determined from the 
fit to the DFT results, are shown as a function of Mn-Mn 
distance in Fig. [IT] The strongest terms are antiferro- 
magnetic terms for exactly those Mn-Mn pairs that (1) 
are approximately 3.5 A to 3.6 A apart, (2) share one 
close O atom bonded to each, and (3) have a Mn-O-Mn 
angle for these bonds less than 125°. These Mn-O-Mn 
pairs involve two Mn(3), Mn(4), or Mn(5) atoms and 
have very unequal Mn-0 distances. Within the Mn-Mn 
pairs the magnitude of the antiferromagnetic interaction 
strictly decreases with increasing length of the longer Mn- 
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FIG. 12: Linkages between Mn in a-Mn2 03 limited to a sub- 
set of those with strong antiferromagnetic interactions as de- 
termined in this work gives the lowest energy collinear mag- 
netic structure on the Mn(3), Mn(4), and Mn(5) sites. 



T .A 




FIG. 13: Similar to Fig. [12] with the Mn-Mn pairs with the 
five strongest AFM interactions shown. One additional in- 
teraction, not shown in Fig. |12| is indicated by dotted lines 
here, and leads to frustrated triangles, which presumably is 
the origin of the noncollinear ground state magnetism. 



O distance. 

Two pair interactions of approximate magnitude 0.085 
eV per pair are nearly degenerate in energy and Mn- 
Mn distance, such that they can not be distinguished 



in Fig. 11 If the Mn-Mn links corresponding to the 3 
strongest terms and one of the 0.085 eV terms are drawn 



(Fig. 12), then the magnetic structure for the Mn(3), 



Mn(4), and Mn(5) sites based on these interactions, all 
antiferromagnetic, is completely determined to be that 



of the computational and experimental collinear ground 
state. On the other hand, if links corresponding to both 



0.085 eV terms are drawn (Fig. 13 1, then there are frus- 



trated triangles where antiferromagnetic interactions can 
not be satisfied for all bonds. A triangle of vector spins 
with frustrated antiferromagnetic interactions has a (pos- 
sible degenerate) ground state with a noncollinear spin 
arrangement; we believe this is the origin of the non- 
collinear antiferromagnetism of a-Mn 2 03. 

We only looked at the ground state magnetic structure, 
but the finite temperature magnetic correlations could 
also be investigated using the DFT model, which would 
allow the nature of possible AFM- AFM transitions in a- 
Mn 2 C>3 (O to be determined. 



Conclusions 

First principles density functional theory DFT+U and 
cluster expansion model calculations, along with inde- 
pendent experimental neutron diffraction structure anal- 
yses, were used to determine the low-temperature crystal- 
lographic and magnetic structure of bixbyite a-Mn2 03. 
Both approaches independently gave nearly identical 
crystallographic and magnetic structures, with identical 
antiferromagnetic ordering along a principal magnetic 
axis and secondary ordering along a single orthogonal 
axis, differing only by a phase factor in the modulation 
patterns. The agreement between the two approaches 
suggests that the ground magnetic state of a-Mn203 is 
largely solved. 

The computational methods exploited a Bayesian ap- 
proach that allows the number of parameters in the clus- 
ter expansion model to exceed the number of input struc- 
ture energies without sacrificing energy predictability. 
The individual magnetic coupling parameters were de- 
termined, showing that specific frustrated antiferromag- 
netic interactions determine the magnetic structure. The 
experimental approach benefited from optimized sample 
synthesis, which produced crystallite sizes large enough 
to yield a clear splitting of peaks in the neutron powder 
diffraction patterns, thereby enabling magnetic-structure 
refinements under the correct orthorhombic symmetry. 

The approaches used here should prove suitable for 
similar problems of magnetic ordering in other complex 
oxides whose magnetic states are determined by a large 
and frustrated set of antiferromagnetic interactions. 
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